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Abstract 

We present a theoretical study of electromagnetic surface waves localized at an 
interface separating a conventional uniform medium and a semi-infinite 1-D photonic 
crystal made of alternate left-handed metamaterial and right-handed material which 
we refer to as left-handed photonic crystal. We find novel type of surface mode's 
structure, the so-called surface Tamm states and demonstrate that the presence of 
metamaterial in the photonic crystal structure allows for a flexible control of the 
dispersion properties of surface states, and can support the Tamm states with a 
backward energy flow and a vortex-like structure. 
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1 Introduction 



The left-handed metamaterial (LHM) with negative permittivity and negative 
permeability has attracted much attention (UH; 0) and triggered the debates 
on the application of the left-handed slab as so-called "superlenses" (0; 0). 
Over 30 years ago, Veselago (6), first proposed that this peculiar medium pos- 
sesses a negative refractive index,which has been demonstrated at microwave 
frequencies in recent experiment (J3j). In such media, there are many interesting 
properties such as the reversal of both Doppler effect and Chernekov radia- 
tion (0), amplification of evanescent waves pj), and unusual photon tuneling 
(0), negative giant Goos-Hanchen effect @). All these phenomena are rooted 
in the fact that the phase velocity of light wave in the LHM is opposite to 
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the velocity of energy flow, that is, the Poynting vector and wavevector are 
antiparallel so that the wavevector, the electric field, and the magnetic field 
form a left-handed (LH) system. 



Interfaces between different physical media can support a special type of local- 
ized waves as Surface waves or surface modes where the wave vector becomes 
complex causing the wave to exponentially decay away from the surface. Aside 
from their intrinsic interest, surface electromagnetics waves have recently been 
proposed (0; H^ l as a way to efficiently inject light into a photonic crystal 
waveguide, or to extract a focused beam from a channel. In periodic systems, 



staggered surface modes are often referred to as Tamm states (lilt |12J) first 



identified as localized electronic states at the edge of a truncated periodic 
potential. An optical analog of linear Tamm states has been described theo- 
retically and demonstrated experimentally for an interface separating periodic 
and homogeneous dielectric media ([lit 14 ) . In optics, the periodic structures 
have to be manufactured artificially in order to manipulate dispersion prop- 
erties of light in a similar way as the properties of electrons are controlled in 
crystals. 



One-dimensional photonic crystal consisting of alternate LHM and conven- 
tional right-handed material (RHM) layers have already been investigated 
through calculating the transmittance or the reflectance the structures (0; 
El) and demonstrating of complete band gap (16). In this letter we demon- 
strate another example of unusual properties of LHM and study electromag- 
netic surface TE waves guided by an interface between RHM and a ID periodic 
structure consisting of alternate LHM and RHM layers which we refer to left- 
handed photonic crystal. We demonstrate that the presence of a LHM allows 
for a flexible control of the dispersion properties of surface states. 



We pay special attention to the terminating layer (or cap layer) of the periodic 
structure and assume that it has the width different from the width of other 
layers of the structure. We study the effect of the width of this termination 
layer on surface states, and explore a possibility to control the dispersion 
properties of surface waves by adjusting termination layer thickness. We find 
novel types of surface Tamm states at the interface with the LH photonic 
crystal. In addition to absolutely novel structure due to sharp jumps at the 
interface of layers resulting from opposite signs of permeability of adjacent 
layers, these modes have a backward flow and a vortex-like structure. The 
surface modes in this case we call left-handed Tamm states. 
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Fig. 1. Geometry of problem. In our calculations we take the following values: 
d\ = 2.5cm, d 2 = 1cm, eq = 1, fio = 1, e± = —4, n\ = —1, e 2 = 2.25, \i 2 = 1 

2 Transfer matrix method and surface waves 



We use the transfer matrix method to describe surface Tamm states that form 
at the interface between an uniform medium and a semi-infinite one dimen- 
sional photonic crystal containing negative refraction metamaterial. Geometry 
of our problem is sketched in Fig. 1. We choose a coordinate system in which 
the layers have normal vector along OZ. As shown in Fig. 1 the crystal is 
capped by a layer of the same material but different width, d c . For conve- 
nience of presentation, we split this termination layer of index of refraction n x 
and width d c into sublayers, of lengths d s + d t = d c . The first sublayers extends 
from z = —d s to z — 0. Then the periodic array that forms the one dimensional 
photonic crystal consists of cells each made of three uniform layers of widths 
d t , d 2 , and d\ — d t whose respective indices of refraction are n 1; n 2 , and We 
consider the propagation of TE-polarized waves described by one component 
of the electric field, E=E y , and governed by a scalar Helmholtz-type equation. 
We look for stationary solutions propagating along the interface with the char- 
acteristic dependence, ~ exp[-iuj(t — (3x/c)), where ui is an angular frequency, 
f3 is the normalized wavenumber component along the interface, and c is the 
speed of light, and present this equation in the form, 
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where k1=u 2 f3 2 / c 2 , both permittivity e(z) and permeability ^l(z) characterize 
the transverse structure of the media. Matching the solutions and derivatives 
at the two interfaces, one finds the elements of transfer matrix, M. 



Ak\d\ 



cosk 2 d 2 + 7- (7— + — ■ ] sink 2 d 2 
2 V k 2 k x , 



M- 



12 



D ik\{d\— 2dt) 



1 «1 . k 2 \ 

- — H smk 2 d 2 

2 \ k 2 k x ) 



(3) 



3 



M 



21 



— iki(di— 2dt) 



- — H sink 2 d 2 

2 \k 2 ml 



M< 



22 



-ik\d\ 



cosk 2 d 2 — - \ -r- H — - | sinkodo 
2 U 2 «J 



(4) 
(5) 



where fej = kynf — ft 2 , i = 1,2 and Ki = ^fei. The parameters /ii and fi 2 are 
permeability of layers in photonic crystal. 

Surface modes correspond to localized solutions with the field E decaying 
from the interface in both the directions. In the left-side homogeneous medium 
(z < —d s , see Fig. 1), the fields are decaying provided ft > EqUq . In the right 
side periodic structure, the waves are the Bloch modes 

E(z) = V(z)exp(iK b z) (6) 



where Kb is the Bloch wavenumber, and ^(z) is the Bloch function which is 
peri odic with the period of the photonic structure (see details, e.g., in Ref. 
(J17|)). In the periodic structure the waves will be decaying provided K b is 
complex; and this condition defines the spectral gaps of an infinite photonic 
crystal. For the calculation of the Bloch modes, we use the well-known transfer 
matrix method dlj 



To find the Tamm states, we take solutions of Eq. (1) in a homogeneous 
medium and the Bloch modes in the periodic structure and satisfy the con- 
ditions of continuity of the tangential components of the electric and mag- 
netic fields at the interface between homogeneous medium and periodic struc- 
ture 0). 



3 Results and discussion 



We summarize the dispersion properties of the Tamm states in the first and 
second spectral gaps on the plane of the free-space wavenumber k = u/c versus 
the propagation constant ft (see Fig. 2) for different values of the cap layer 
thickness d c . 

As mentioned above, Tamm states exist in the gaps of the photonic bangap 
spectrum (unshaded regions in Fig. 2). An important point that we can see 
from Fig. 2 is that the sign of the slope of the dispersion curve changes at 
one (or two) extremum point at first (or second) bandgap, respectively. The 
slope of the dispersion curve determine the corresponding group velocity of the 
mode. According to the energy flow of modes (see Fig. 4) positive (negative) 
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Fig. 2. Dispertion properties of the Tamm states in the first and second spectral 
gaps. Shaded: bands of one-dimensional LH photonic crystal. Solid: dispersion of 
the Tamm states with LHM. Corresponding values of the cap layer thickness are 
indicated next to the curves, Lower and upper insets show blow-up regions of small 
j3 for d c = 0.6cfi and d c = 0.35di respectively. Points (a), (b), and (d) correspond 
to the mode profiles presented in Figs. 3(a, b, d). 

slope of dispersion curve, here, correspond with backward (forward) energy of 
mode which we call backward (forward) Tamm states. So, for a thin cap, d c = 
0.6c?!, the backward Tamm states occur at the incidence angles corresponding 
(3 > 1.17, while backward Tamm state for thick cap d c = 1.2d\ take place at 
wide ranges (3 > 1.02. Meanwhile the LH Tamm surface wave has a vortex-like 
energy flow pattern. 

As a result of different slope of the dispersion curve of the LH Tamm states, we 
observe the mode degeneracy, i.e, for the same frequency u (or wavenumber 
k), there exist two modes with different value of f3 in first bandgap (see lower 
inset in Fig. 2) while, we can see three modes with different value of (3 in 
second bandgap (see upper inset in Fig. 2). There are modes with negative and 
positive group velocities (with respect to the propagation wavevector). Such 
modes are termed as forward and backward , respectively. In the forward wave, 
the direction of the total energy flow coincides with the propagation direction, 
while in the backward wave the energy flow is backwards with respect to 
the wavevector. Physically, this difference can be explained by looking at the 
transverse structure of modes. In Figs. 3 (a,b) we plot the profiles of the 
two modes having the same frequency k = uj/c = 1.037cm -1 , with different 
longitudinal wavenumber (3. Corresponding points are shown in the lower inset 
in Fig. 2. For the mode (a), the energy flow in the homogeneous RHM exceeds 
that effective backward energy flow in the periodic structure (slow decay of 
field for z < —d s and fast decay into the periodic structure), thus the total 
energy flow is forward. For the mode (b) we have the opposite case, and the 
mode is backward. To demonstrate this, in Fig. 4 we plot the total energy 
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Fig. 3. Examples of the LH surface Tamm states, (a) Forward LH mode, d c = 0.6<ii, 
k = 1.037cm- 1 , = 1.05; (b) Backward LH mode, d c = 0.6di, k = 1.037cm" 1 , 
P = 1.36; (c) Guided backward LH mode, d c = 3di, k = 0.9965cm" 1 , p = 1.2; (d) 
Backward LH second bandgap mode, d c = 0.35di, k = 2.029cm" 1 , p = 1.?; 
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Fig. 4. Total energy flow in first (solid) and second (dashed) bandgap Tamm modes 
vs. p. Corresponding values of the cap layer thickness are indicated next to the 
curves. 

flow in the modes as a function of the wavenumber P for first and second 
bandgap surface modes. These results confirm our discussion based on the 
analysis of the dispersion characteristics. One can see from Fig. 4 that the 
increasing the cap layer thickness, d c leads to wide range of backward modes 
(campare energy flow of d c = O.Qdi and d c = 1.2g?i). In the conventional 
geometry consisting of RHM layers hence right-handed (RH)Tamm state the 
energy flow in all parts of the wave are directed along the wave number and 
the localized surface waves are always forward. An example of the second-gap 
mode profile is shown in Fig. 3(d). One can see that the second band modes are 
generally weaker localized at the interface then the modes from the first band 
gap. In addition, the oscillation period of the transverse structure of modes at 
the first bandgap is twice of the physical structure of photonic crystal while 
these are equal at the second bandgap. 



Finally, in Fig. 5 we plot the existence regions for the surface Tamm modes for 
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Fig. 5. Existence regions of the surface Tamm modes for first bandgap (a) and 
second bandgap (b). The modes do not exist in the shaded regions. 

the both first (Fig. 5(a)) and second (Fig. 5(b)) bandgaps on the parameter 
plane (d c , 0). Shaded is the area where the surface modes do not exist. Increas- 
ing the cap layer thickness, d c , we effectively obtain a dielectric waveguide, one 
cladding of which is a homogeneous medium, while the other one is a LH pho- 
tonic crystal. In such typical mode is shown in Fig. 3(c). Fig. 5 shows 
that the non-existence regions of Tamm states qualitatively different for first 
and second bandgaps. While there exist three discrate non-existence regions 
of Tamm states for first bandgap, we found five discrate non-existence regions 
of Tamm states for second bandgap including backward and forward modes. 



4 Conclusions 



We have studied novel types of electromagnetic surface waves guided by an 
interface between a homogeneous conventional material and a one-dimensional 
LH photonic crystal consisting of alternate LHM and RHM layers. We have 
shown that in the presence of a LHM in the photonic crystal structure the 
surface Tamm waves firstly have a absolutely novel shape with a sharp jumps 
in the interface of layers due to opposite signs of the permeability in two adja- 
cent media, secondly can be either forward or backward. It must be noted that 
Tamm states for conventional structure are always forward. We have studied 
the properties of the first and second bandgap Tamm states and analysed the 
existence regions for both band modes. We believe our results will be useful 
for a deeper understanding of the properties of surface waves in plasmonic and 
metamaterial systems. 
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